Hyperbolicity and the effective dimension of spatially-extended dissipative systems 
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Using covariant Lyapunov vectors, we reveal a split of the tangent space of standard models of 
one-dimensional dissipative spatiotemporal chaos: a finite extensive set of A'^ dynamically entan- 
gled vectors with frequent common tangencies describes all the physically relevant dynamics and 
is hyperbolically separated from possibly infinitely many isolated modes representing trivial, expo- 
nentially decaying perturbations. We argue that N can be interpreted as the number of effective 
degrees of freedom, which has to be taken into account in numerical integration and control issues. 
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Nonlinear dissipative partial differential equations 
(PDEs) are ubiquitous in the description of pattern- 
forming, chaotic, and turbulent systems [1]. Even though 
they are formally infinite-dimensional dynamical sys- 
tems, it is now well accepted that their chaotic solutions 
evolve in an effective manifold of finite dimension. For 
many generic PDEs such as the Kuramoto-Sivashinsky 
(KS) and the complex Ginzburg-Landau (CGL), it is in 
fact proven that trajectories are first exponentially at- 
tracted to a finite-dimensional invariant manifold called 
the inertial manifold |2||. This object, however, remains 
largely formal, as there does not exist a constructive way 
of determining of which modes it is composed. Similarly, 
trajectories eventually fall into a global attractor of fi- 
nite Hausdorff dimension. For large systems, this dimen- 
sion, as well as other quantities measuring the amount 
of chaos in the system, can be estimated via the calcula- 
tion of Lyapunov exponents. These dimensions remain, 
however, global quantifiers. One approach to determine 
which modes actually compose and contribute to the dy- 
namics is that pursued, e.g., by Cvitanovic et al. [3|, but 
it is difficult and limited to rather small systems. 

A related difficulty lies in the numerical integration 
of dissipative PDEs (which remains the primary way 
of studying their often chaotic solutions). If the finite- 
dimensionality of their attractors justifies that a numer- 
ical study is possible at all, there is no a priori criterion 
to define the minimal resolution for a faithful simulation, 
and in practice, one typically checks the convergence of 
results upon increasing the resolution. 

In this Letter, we use covariant Lyapunov vectors 
(CLVs), recently made numerically accessible thanks to 
an efficient algorithm [J], to show that the tangent dy- 
namics of large KS and CGL systems is essentially char- 
acterized by a well-defined set of "physical" modes. Be- 
cause the covariant vectors span the intrinsic (Oseledec) 
subspaces corresponding to each Lyapunov exponent, 
and thus allow access to hyperbolicity properties, we are 
able to show that the physical modes are decoupled from 



the remaining set of hyperbolically "isolated" degrees of 
freedom. In the context of dissipative partial differential 
equations, our results imply that a faithful numerical in- 
tegration needs to incorporate at least as many degrees of 
freedom as the number of such physical modes and that 
further increasing the resolution increases the number of 
degrees of freedom associated with the trivially decaying 
isolated modes. 

We first focus on the one-dimensional KS equation, 
taken here as a prototypical dissipative PDF showing 
space-time chaos [l|, l5|. It governs a real field u{x,t) 
according to 



dtu 



-d^u — d^u ~ udxU, X e [0, L\. 



(1) 



Figure [T] shows the Lyapunov spectrum for a fixed sys- 
tem size L = 96 but different spatial resolutions and pe- 
riodic or rigid boundary conditions (PBC or RBC) [6]. 
The spectrum consists of two parts: first a smooth re- 
gion of positive, zero, and some negative exponents, then 
a rather steep region of negative exponents arranged in 
steps of two for PBC. The two regions are separated by an 
abrupt change in slope (bottom inset). Remarkably, the 
spectra for different spatial resolutions overlap with the 
extra exponents coming from the higher resolution simply 
accumulating in the second region, at the negative end 
of the spectrum (upward and downward triangles). Thus 
the threshold index separating the two regions stays un- 
changed (here around j — 40) upon increasing resolution. 
Note also that the boundary conditions only change the 
multiplicity of modes in the second region, where every 
other mode is exactly the same for both PBC and RBC. 
The above observations suggest that modes in the sec- 
ond region, hereafter called "isolated" , for reasons given 
below, are residual, highly damped degrees of freedom 
not necessary to describe properly the essential dynam- 
ics. In contrast, the modes in the first region, which 
should be intimately associated with phase space dynam- 
ics, will be called "physical" . In the following, we sub- 
stantiate this intuition on a rigorous basis studying the 
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FIG. 1: (color online) Lyapunov spectrum A'--'' made of ex- 
ponents arranged by decreasing value for the one-dimensional 
KS equation with L — 96. Black upward triangles: fccut = 
42-2n/L and PBC. Red downward triangles: fccut = 85-2tt/L 
and PBC. Blue dots: fccut = 170 ■ tv/L and RBC {u{L,t) = 
u{0,t) = 0). In all three cases the spectrum yields a positive 
region with maximum Lyapunov exponent A'^' ~ 0.09. Top 
inset: A"' vs j* for the large-j (isolated) modes. Bottom 
inset: close-up around the threshold. 



CLVs associated with Lyapunov exponents. The CLVs 
for the isolated modes, contrary to those of the physi- 
cal modes, possess an approximately sinusoidal, delocal- 
ized structure as can be seen from their power spectra 
(Fig. [^b) , as well as directly from their rather uneventful 
spatiotemporal evolution (Fig. [2ja,). With PBC, the two 
modes forming a step show the same dominant wavenum- 
ber with an arbitrary phase shift. This is not the case of 
RBC where the phase of the sinusoidal structure is fixed 
at the boundary. The peak wavenumber k^J^^, which lin- 
early increases with j (top panel of Fig.[2j:), is in fact just 
the j-th wavenumber allowed in the given spatial geome- 
try (multiplicity taken into account): fcpj,g^]^ = [j/2]-27r/L. 

For large-enough j, A^^'' ~ "(^icak)'' '^ ^J^ i^'^P inset 
of Fig. [1]), which indicates that the values of the Lya- 
punov exponents of the isolated modes are governed by 
the stabilizing linear term of the KS equation (i.e. the 
fourth-order derivative). 

The sinusoidal structure of the isolated modes indi- 
cates that they are nearly orthogonal to each other. In- 
deed, distributions of the angle 6 between pairs of CLVs 
of indices j > 42 are peaked at 7r/2 and drop rapidly 
near and tt. This is also true if only one of the vectors 
is taken in the region j < 41, but for any pair of vec- 
tors taken from this region, the angle distribution spans 
the whole [0,7r] interval (Fig. [2ji). A careful analysis of 
the angle distributions reveals that those involving iso- 
lated modes seem to have an essential singularity near 
(and tt): p{e) - exp(-const./6') (Fig. H?). Given the 
sharpness of this behavior, we cannot conclude about the 
possibility for these distributions to be strictly bounded 
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FIG. 2: (color online). Properties of CLVs for the KS system 
(L = 96, fccut = 42-27r/L, and PBC). (a) Spatiotemporal plot 
of a typical vector in the isolated region {j = 46, total time 
100). (b) Spatial power spectra of vectors of indices j = 1, 
16, 32, 38, 44, 52, 60, 68, 76, 84 (from left to right in peak 
position), (c) Top panel: peak wavenumber in the power 
spectra (red circles) and fc = [j/2] ■ 2n/L (black line) . Bottom 
panel: DOS violation fraction j/.f"''"' for pairs of neighboring 
vectors (pairs within the same step are omitted), (d) Angle 
distributions between pairs of vectors, (e) Same as (d) but 
with a different abscissa. 



away from zero. The above results determine accurately 
the threshold at j = 41, and therefore the number of 
physical modes, and give the precise definition of physi- 
cal and isolated modes: contrary to the physical modes, 
the isolated modes do not have any tangencies with other 
isolated modes (except the partner in the same step) nor 
with physical modes. They can be said to be "hyperbol- 
ically isolated." 

The absence of tangency for isolated modes can be con- 
firmed from another viewpoint, that of the so-called dom- 
ination of the Oseledec splitting (DOS) [Ti], which quan- 
tifies, loosely speaking, the degree of dynamical isolation 
of the Oseledec subspaces from each other due to the 
strict ordering of Lyapunov exponents. Let Xr (t) be the 
finite-time Lyapunov exponent averaged over a period r 
around time t. The splitting of the space formed by the 
vectors associated with the modes ji and J2{> ji) is said 
to be dominated if Xi^'\t) > xi^^\t) holds for all t with 
T larger than some finite To. It is mathematically proven 
that DOS implies absence of tangency between the Os- 
eledec subspaces, or the CLVs [7]. To quantify DOS, we 
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define, following [| AX^.^''^'\t) = 
with jmax = niax(ji,J2) and j„ 
measure the time fraction of DOS violation Vr'''''^ = 
(0(AAr^'"'^ (t))}, where Q{z) is the step function and 
(• • • ) denotes the time average. The result is shown in 
the bottom panel of Fig. [2j; for pairs of neighboring ex- 
ponents: i^T '"* with r ~ 0.2 drops sharply near the 
threshold and becomes strictly zero for j > 41. In fact, 
z/V'"' stays zero for any pairs with j or j' > 41, though in 
some cases slightly larger values of r are required. This 
confirms the absence of tangencies of the isolated modes, 
already seen from the angle distributions, and the num- 
ber of physical modes. 

We now turn our attention to a different case in or- 
der to test the validity of our results beyond the simple 
KS equation. Let us consider the CGL equation, whose 
universal relevance and genericity is now well-established 
[l|, |9|. In one space dimension, it governs a complex field 
W{x,t) according to: 

dtW = W - {I + iP)\W\^W + {1 + ia)d'^W . (2) 

In the following we consider a so-called "amplitude tur- 
bulence" regime [10|], i.e. a strongly chaotic regime 
where amplitude and phase modes evolve on rather short 
time- and length-scales. (Results for other regimes, such 
as phase turbulence, will be presented elsewhere [llj.) 
Specifically, we use a = —2.0, /? = 3.0, L = 64, and PBC. 

For sufficiently high spatial resolution, the Lyapunov 
spectrum indeed shows an isolated, stepwise region as 
for the KS equation (Fig. [S]:), but here the multiplic- 
ity of each step is four and the spatiotemporal evolution 
of isolated modes reveal patches of traveling waves (Fig. 
[3^b). (The effects of spatial resolution and/or bound- 
ary conditions are also similar to our observations on the 
KS equation.) As for the KS equation, this is in agree- 
ment with the linear stability analysis: normal modes 
are traveling waves of the form exjp[i{±kx — ujkt) + At] 
with A = I — k^ and ujk = ak^. Indeed, the values of 
the Lyapunov exponents in the isolated region behave 
like this and the vectors are composed of traveling waves 
propagating at a velocity of ±tj/fc. Compared to the KS 
case, the additional multiplicity of two comes from the 
degeneracy between k and —k modes. Moreover, this de- 
generacy implies that these two modes are in fact mixed 
up in a single vector: isolated vectors are either in the 
pure k mode (Fig. [3^), in the pure —k mode, or patches 
of the two (Fig. [3)). 

In spite of these differences, the isolated modes for 
CGL remain dynamically isolated from any other mode. 
Although Vr''' measured with small r is not zero in 
the isolated region, a clear threshold is found when con- 
sidering larger r values: Contrary to what happens for 
ji,J2 < 86, all I'r^'''^' for J2 > 87 show a decrease faster 
than exponential, indicating the existence of a finite r 
beyond which Vr^'-''^' is zero (Fig. |3Ji). Consequently, we 
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FIG. 3: (color online) CGL equation in the amplitude tur- 
bulence regime (L = 64, fccut = 31 ■ 2n/L, and PBC). (a,b) 
Spatiotemporal plots of the phase component of a typical vec- 
tor j = 91 in the isolated region (same trajectory, but at two 
distant periods of time, during a total time of 20 for each 
plot), (c) Lyapunov spectrum; inset: close-up around thresh- 
old, (d) Time fraction z^t '"' of DOS violation, as a function 
of T {j = 78, 82, 86, 90, 94, from top to bottom). 



can also define the isolated modes by their DOS as for 
the KS equation, and determine exactly the number of 
physical modes (here at j = 86). Indeed, physical modes 
have no tangencies with any isolated modes {j > 87). 
Moreover, angle distributions involving isolated modes 
show, like for the KS equation, the essential singularity 
^ exp(— const. /0) near tangency. 

We now discuss our results. We first note that the 
threshold separating physical from isolated modes does 
not coincide with the appearance of steps in the Lya- 
punov spectrum (for PBC). Lyapunov exponents alone 
can only provide a good guess: For the KS and CGL sys- 
tems treated above, the first steps start respectively at 
j — 40 and j = 79, whereas the exact thresholds are at 
j = 41 and j = 86. Indeed, a closer scrutiny of the expo- 
nents reveal that the first steps are actually not perfect. 

Let us now specify the implication of the lack of tan- 
gencies for the isolated modes. Suppose that we add 
to the dynamics an infinitesimal perturbation along the 
CLVs of some isolated modes. Then this perturbation 
decays exponentially to zero as indicated by their nega- 
tive Lyapunov exponents, and the absence of tangencies 
implies that this does not induce any perturbation along 
directions spanned by the other Lyapunov modes. In 
contrast, perturbations along physical modes will propa- 
gate to other physical modes through tangencies between 
them, and could eventually induce activity in the modes 
associated with positive exponents, growing to consid- 
erably affect phase space dynamics even if the initial 
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FIG. 4: (color online). Extensivity of the Lyapunov spectrum 
for the KS equation with PBC. Inset: quantities indicating ef- 
fective dimensions of the system: number of non-negative ex- 
ponents (circles), Kaplan- Yorke dimension (squares), metric 
entropy (diamonds, multiplied by 50), and number of physical 
modes (triangles). 



perturbation was made in the direction associated with 
negative exponents. In this sense, the dynamics corre- 
sponding to the physical modes is highly entangled, but 
completely decoupled from the decaying dynamics of the 
isolated modes. Therefore, all the degrees of freedom 
associated with physical modes are necessary to faith- 
fully describe phase space dynamics, while adding further 
more degrees of freedom associated with isolated modes 
does not affect phase space dynamics in any significant 
way. The number of physical modes scales linearly with 
the system size L: rescaled Lyapunov spectra collapse 
both in the physical and in the isolated region with the 
stepwise structure retained (Fig. 2]). In particular, the 
dimension density defined from the number of physical 
modes is larger than others, for instance being almost 
twice the Kaplan- Yorke dimension [iJl (inset of Fig. [4]). 
It is therefore natural to interpret the number of physical 
modes as an embedding dimension of the global attrac- 
tor. Furthermore, we speculate, pending mathematical 
rigor, that the number of physical modes could be re- 
lated to the dimension of the inertial manifold, which is a 
positively invariant and exponentially attracting smooth 
manifold embedding the global attractor (the subspace 
spanned by the physical modes would then be the lo- 
cal linear approximation of the inertial manifold [13j). 
While this conjecture would drastically lower existing es- 
timates for the dimension of the KS inertial manifold 
D < const. X L^-^^ [l3|, it is however in line with the 
extensivity of chaos already observed in the past for the 



same models |l4ll5|. 



Our results suggest that a faithful numerical integra- 
tion of PDEs needs to incorporate at least the degrees 
of freedom associated with all physical modes. More- 
over, the extensivity of the number of physical modes im- 
plies that the minimal resolution evaluated in systems of 
moderate size carries over to arbitrarily large system size 
L. Our obtained number of effective degrees of freedom 
would also be helpful to determine the minimal number 
of constraints necessary for a full control of a continuum 
system, with applications to real situations like e.g. the 



suppression of ventricular fibrillation. 

In summary, we have shown, using Lyapunov analysis, 
that the tangent space of two representative nonlinear 
dissipativc PDEs systems can be divided into two parts: 
a finite dimensional manifold spanned by strongly inter- 
acting physical modes, and the remaining set of isolated, 
strongly damped modes. We demonstrated that isolated 
modes are hyperbolically separated from all other modes, 
and thus satisfy the property of domination of Oseledec 
splitting. Similar results were obtained also for a chain of 
diffusively-coupled tent maps (not shown, 11]). We have 
interpreted the number of physical modes as an embed- 
ding dimension of the global attractor. The extensivity 
of this dimension could also be of interest in view of the 
studies by Egolf et al. [l5| arguing about the "building 
blocks" of spatiotemporal chaos. We hope our results will 
trigger work to clarify these issues at the mathematical 
level, as much as numerical investigations of other dissi- 
pativc systems, like those yielding fully developed turbu- 
lence, for which no rigorous proofs are known about the 
existence of an inertial manifold. Besides their theoreti- 
cal importance, our results are also useful for the numer- 
ical integration and control issues of dissipativc PDEs. 
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